Ph! 



X 



Jan Lorenz, Stefano Battiston: 
Systemic risk in a network fragility model analyzed with 
probability density evolution of persistent random walks 



Systemic risk in a network fragility model analyzed with 
probability density evolution of persistent random walks 

r^ I Jan Lorenz, Stefano Battiston 

Q ■ Chair of Systems Design, ETH Zurich, Kreuzplatz 5, 8032 Zurich, Switzerland 
CN . Ij alorenzQethz . ctil IsbattistonOethz . chl 

>' 



Abstract 



o 

Ps| I We study the mean field approximation of a recent model of cascades on networks rel- 

evant to the investigation of systemic risk control in financial networks. In the model, the 
hypothesis of a trend reinforcement in the stochastic process describing the fragility of the 

r^ [ nodes, induces a trade-off in the systemic risk with respect to the density of the network. 

Increasing the average link density, the network is first less exposed to systemic risk, while 



O I above an intermediate value the systemic risk increases. This result offers a simple explana- 

^ ' tion for the emergence of instabilities in financial systems that get increasingly interwoven. In 

lyj . this paper, we study the dynamics of the probability density function of the average fragility. 

O I This converges to a unique stable distribution which can be computed numerically and can 

c/3 ■ be used to estimate the systemic risk as a function of the parameters of the model. 

1 Introduction 

> 

t;;-|- ■ 1.1 Systemic Risk in Financial Networks 

oo ■ 

'^ ' A network of interdependent units which, individually, are susceptible to fail, is potentially 

exposed to multiple joint failures of a significant fraction of units in the system. This is the notion 
that is usually associated with the term systemic risk. Systemic risk is particularly important 

^ ■ in the context of infrastructure networks, such as power grids, and in financial networks. These 

latter should be meant in a broad sense, including units of different types, such as business firms, 
insurance companies, banks, mutual funds and other financial institutions that are linked by 
credit relationships. For instance, if one or more firms fail and are not able to pay back their 

5^ I debts to the bank, this affect the balance sheet of the bank which might try to improve its own 

situation by increasing the interest rate to the other firms, causing other failures among the 
firms. If finally the bank itself fails, this affects negatively the banks that are linked to it by 
interbank loans. This is somehow similar to failure cascades in power grids where a failing power 
line implies a higher load an other lines which might bring them to fail. The size distribution of 
such failure avalanches is one way of quantifying the systemic risk. 

There is a growing body of literature in economics on financial networks, that investigates also 
the issue of systemic risk. While banks-firms credit relationships have been extensively studied 



IM 



Jan Lorenz, Stefano Battiston: 
Systemic risk in a network fragility model analyzed with 
probability density evolution of persistent random walks 



(for an overview, see p. 7]). only recent works have analysed phenomena of financial contagion in 
interbank credit [21 [12] and trade credit. The latter, is a form of credit among business firms, 
typically in a supplier-customer relation, which has been less investigated despite the fact that in 
some countries it represents a significant part of the short-term liabilities of the corporate sector 
[9]. In the literature on complex networks only few works have dealt with financial networks, 
mainly in the context of self-organized criticality [HIH]. Most of those works suggest that when 
the degree of the nodes in the network increases the network is less exposed to systemic risk. 
In some cases, the evidence that systemic failures may more rare but also more severe has been 
found (see for instance |14|). 

1.2 The Fragility Model for Cascades on Networks 

In this paper, we consider the model of cascades on networks introduced by [6], in which a clear 
tradeoff emerges in the systemic risk, as a function of the network density. This means that up to 
an intermediate level of network density there is a benefit in creating links between units because 
they allow to diversify the risk. However, above a certain level of density, the probability of many 
joint failures increases. This effect depends on the presence of a sort of trend reinforcing term 
in the dynamics of the fragility of the nodes. The fragility is a state variable that determines 
the failure of the node, when it exceeds a given threshold, as well as subsequent transfer of 
some damage to the connected nodes. The trend reinforcing of the fragility corresponds to the 
following idea. If the fragility of a firm at the end of the year has reduced compared to last year, 
the firm is rated better in terms of solvency and it has easier access to credit. Conversely, if the 
fragility has increased, the firm faces worse conditions for credits and thus additional cost that 
are likely to increase its fragility furthermore. Notice that, through the links in the network, 
this propagates also to the neighbours, since the fragility of the firm affects the fragility of the 
neighbours. For instance, hedge funds leverage even small differences in performance across firms 
by 'short-selling' the stocks of the slightly worse ones and 'going long' on the slightly better ones. 
Thus, even small differences in the evolution of two firms may matter a lot. Further on, effects 
like predatory trading [TO] may induce trend reinforcing. 



1.3 Outline of the Paper 

In [^, some analytical results supporting the simulations are found, based on separating the pro- 
cess of the evolution of fragility (approximated as a time-dependent Ornstein-Uhlenbeck process) 
and the cascade process (where the size of an avalanche is expressed as the fix point of an equa- 
tion for the number of failures). Here, we provide an alternative analysis of the tradeoff regarding 
systemic risk mentioned above. We consider the stochastic process defined by the mean field ap- 
proximation of the fragility of the individual node. This is now a stochastic process for a single 
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variable, and it is also clear that, having reduced the system to one single variable, the cascading 
part of the process is excluded by construction. In this approximation the failure probability can 
be taken as a proxy for systemic risk. In fact, the mean field approximation is valid when all 
units behave in a similar way. We study some mathematical properties of the process and we 
provide a simple method to show the existence of a tradeoff in systemic risk as function of the 
density of the network. The method is based on recognizing that the process is a combination of 
a Gaussian Random Walk (RW) and a Persistent Random Walk (PRW). PRW ^9\ is a variant 
of the classic RW in which the walker has a probability p to keep the direction of his former 
movement and 1 —p to switch direction. The process is sometimes called correlated random walk. 
It is approximated by the Telegraph's equation [TH1IT3] in the limit of continuous time and space. 
It differs from RW in the scaling with time of the variance of the displacement of the walker. In 
our model, the dynamics in time of the fragility induces a dynamics on the probability density 
function of its values. This dynamics has an exact analytical expression and the systemic risk is 
measured as the number of failures in the stable distribution of fragility. It is possible to proove 
the existence, uniqueness and convergence to a stable distribution, based on the Birkoff-Jentzsch 
theorem which extends the Perron- Frobenius Theorem to infinite dimensional vector spaces. We 
cannot provide an closed-form expression of the sytemic risk as a function of the parameters of 
the model, but we compute the systemic risk numerically, by iterating the dynamics on the pdf. 
We show in this way that the systemic risk has indeed a minimum as function of the network 
density. 

The paper is organized as follows. In Section [2] we introduce the model. In Section [3] we analyze 
the model: first, we describe the mean-field approximation of the dynamics and we show how it 
can be described by using a PRW. Then in Section [321 we derive the dynamics on the probability 
density function and we prove existence and uniqueness of the stable pdf. In Section H] we report 
the results of the numerical computation of systemic risk. In Section O we check the robustness 
of our results with respect to the type of noise that enter in the stochastic process of the fragility 
and some other slight modifications. In Section [6] we summarize the results and we draw some 
conclusions. 



2 The model 

In this section, we describe the network fragility model introduced in [6]- Consider a set of n firms 
connected in a network, each associated with two state variables, the size a and the fragility if. 
The first captures the notion of a proxy for the size of the firm, such as its output. The fragility 
captures the notion of financial fragility of the firm. This is measured for instance in terms of 
its net worth: when the net worth decreases down to zero, the firm is not able to pay back its 
debts and goes bankrupt. So the larger the net worth, the smaller the fragility. As shown in [3], 
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in a network of firms linked by supply-customer relationships, the net worth of a firm evolves 
as a stochastic process that depends on the net worth of the neighboring firms. The interaction 
with the neighbors results in an averaging term and in a trend reinforcing term. Each firm has 
a portfolio of suppliers and customers, which reduces the impact of the fluctuations of prices 
and shocks both from the suppliers and customers, thus resulting in the averaging term. On the 
other hand, if the production cost increases when the net worth of the firm and its neighborhood 
is decreasing (because it is more costly for the firm to access the credit it need for production), 
this results in a trend reinforcing term |5]. Following [6] we model directly the fragility of firms 
as a stochastic process confined in the interval [0,0], where 9 is the failure threshold. 

Firms are connected in a weighted and directed graph with adjacency matrix W E M"^". W is 
non-negative and row-stochastic (i.e. ^ • Wij = 1). 

As a first step, let us look at the following equation for the evolution of the fragility of the set of 
firms 

^{t + 1) = Wifit) = vfV(o) (1) 

where ip = [c^i, . . . , c^„] is the vector of fragility values. If Wij is positive, then the fragility of 
firms j contributes to a fraction Wij to the value in the next time step of the fragility of firm i. 
In other words, the fragility of firm i at time i -|- 1 is a weighted arithmetic mean of the fragility 
values of the neighboring firms. Under some conditions about connectivity in the network, the 
values of fragility of the firms will converge in time to a same value — namely if the matrix W has 
only one essential class of indices which is primitive (this is shown in [16], such matrices are called 
regular if they are row-stochastic, as in our case). If there are more then one essential classes 
the fragilities in these classes converge internally to the same value, as well as all inessential 
firms which have connections exclusively to this essential class. But there is no interplay with 
fragilities in other essential or inessential classes. If an essential class is not primitive there is 
some internal cycling of fragility values. See [HI [7] for the results in the context of conditions 
of finding consensus in a group of experts. So, for graph with high link density we could assume 
that the fragility values will converge to the same value. 

We now introduce additive stochastic shocks and trend reinforcing. 

^{t + 1) = W{ip{t) + am) + asign{W{ip{t) - ip{t - 1))) (2) 

In the equation above ^(t) is a vector of iid random variables, ^i(t), . . . ,^„(i), drawn from a 
distribution /^, with expected value zero and standard deviation one and no skewness (i.e. its 
probability density function is symmetric). The parameter a determines the standard deviation 
of shocks and is also called the noise level. The fragility of each firm receives, as a net shock, the 
weighted average of the shocks that hit the fragility of the firms in its neighborhood. In other 
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words, the firm hedges the risk for upward shocks to its own fragihty, by sharing the shocks 
with other firms. In the second term of the equation, the sign is apphed component-wise (for 
completeness we define 0(— 1) = 0) and a is a constant that we call the trend strength. A fixed 
constant a is added if the difference between the current average fragility in the neighborhood 
and that at the previous time step is positive (i.e. if fragility has increased) and is subtracted if 
the difference is negative (if fragility has decreased) . 

As a result of the dynamics of Eq. ^, the values of fragility may very well go out of the interval 
[0,9]- Therefore, (pi{t + 1) is set to zero if (/>j(i + 1) ^ [0,^] . For firms whose fragility would 
go below zero this means that their fragility cannot become lower than that. For firms that get 
above 9 this means that they go bankrupt and are replaced by a new firm with initial fragility 
zero. So, Eq. ([2]) can be stated as 

^{t + 1) = l[o,e] {W{^{t) + a^t)) + asign{W{^{t) - ^{t - 1)))) 

where l[o,6»] is the (componentwise) indicator function (e.g. Ijo^ejCv') = 1 if v' ^ [Oi^] ^i^d 
otherwise, also known as X[o,6»])- 

In the following we will omit l[o,6»] when we describe dynamics because the reset to zero when a 
firm fails is not the only reasonable choice. We discuss some variations at the end of the paper. 
In any case throughout we assume that the the process is somewhere reset when it gets out of 
[0,0]. 

In the original model in [6], when a firm i goes bankrupt, some damage, proportional to the 
size Qi of the firm is transferred to the fragility of neighbors. If, as result, the fragility of some 
neighbors exceed the threshold 9, they, in turn, transfer a damage to their (surviving) neighbors. 
This cascading process occurs at a faster time scale than the dynamics above. In this paper, we 
do not use at all the cascading part of the model. So Eq. ([2]) describes completely the dynamics 
we study here. 

3 Model analysis 

Since the dynamics depends on the relative magnitude of the parameters a, a and 9. we can fix 
9 = 1 without loss of generality. For abbreviation we define the difference A(/j(i) = ip{t) — ip{t — l). 

If W is the unit matrix (i.e. there is no hedging of risk) ^ reduces to 

^i{t + 1) = ^i{t) + aii{t) + asign(Ac/5i(t)) (3) 

for all i. 
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If firms are connected in a complete graph and share their fragility shock to an equal proportion 
with all other firms, then Wij = ^ for all i,j. In this case, the fragility of each firm, evolves as 
the average (p{t) = ^ X^j=i ^i{t)- Then, the central limit theorem implies that 

^{t + 1) = 0(t) + -^C{t) + asign(A0(t)). 



In general, if each firm is connected, on average, to k < n other firms, one can make a mean-field 
approximation of the dynamics of the fragility of each firm and write 

0(t + 1) = 0(t) + -^e(t) + asign(A0(t)). (4) 

The parameter k is the average number of hedging partners or hedging level. In other words, the 
stochastic process on (j) represents the evolution of the average fragility of the economy where each 
firm has on average k hedging partners. In this approximation, increasing the average number of 
hedging partners k decreases the standard deviation of the shocks a by a factor of \/k. Intuitively, 
one can expect that the failures becomes less frequent, because, the smaller are shocks at each 
time step, the longer it takes to eventually hit the threshold 6. However, if the noise level a is 
very small compared to the trend strength a, the second term in Eq. ([5]) dominates. In particular, 
if the fragility was increasing from time t — 1 to time t, then the second term is for sure equal 
to +a while the first is probably very small and therefore the fragility will also increase at time 
t + 1. Therefore, the noise level or equivalently, the average number of neighbors in the network, 
seems to play a crucial role for the probability of a given firm to hit the fragility threshold. 

As an example. Figure [1] shows six trajectories of the stochastic process defined in Eq. ([5]) for a 
fixed value of trend strength a and decreasing value of noise level a. 

In the following, we will investigate the role of noise on the probability of failure by computing 
the pdf of (j) in the limit of large t, which represents the probability distribution of fragility in the 
steady state of the process. Such pdf can be interpreted both as the firm's individual probability 
of having a given value of fragility and as an histogram of fragility values of an ensemble of firms. 

3.1 Dynamics of Fragility as Persistent Random Walk 

Since varying the hedging level k is equivalent to varying the noise level, in the following we 
definitely drop k from Eq. ^ and we study the process 

^{t + 1) = 0(t) + a^{t) + asign(A</.(t)) (5) 

Assuming that the boundary conditions are not effective during two consecutive time steps, we 
can derive from ([5]) the expression of (j){t + 2) in terms of (j){t). 

0(t + 2) = <p{t) + a(e(t + 1) + m) + « [sign(A0(t)) + sign(a^(t) + asign(A(/)(t)))] . (6) 
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a = 0.05, o = 0.2, 21 bankruptcies 



a = 0.05, = 0.1,11 bankruptcies 




100 



100 



Figure 1: Example of six trajectories of the stochastic process for 4> with fixed trend strength 
a and decreasing noise level a. The number of failures first decreases but then increases. Since 
decreasing the noise level is equivalent to increasing the hedging level, the figure suggests that 
there is an optimal hedging level which minimizes the number of failures. 



Obviously, the last term in the square parentheses can only take the values —2, or 2, depending 
on the sign of A(/)(t) and the probability 



This probability is 



Pr(sign(o-^(t) + asign(A(/)(t))) = sign(A(/)(t))). 



Pr(crC < a) 



fa£_{x)dx = I f^{x)dx 



due to the symmetry of /^. We define q{a, a) := Pr(o"^ < a) as the probability to keep the trend. 
Denoting with F^{x) the cumulative distribution function (cdf) of ^ then it holds q{a, a) = F^{^). 
We can then reformulate the process ^ as 



(pit + 1) = 0(t) + aC{t) + atr(f ) 



(7) 



where (/)(t + 1) is set to zero if it falls out of the interval [0,0]. The function 'tr' is the discrete 
stochastic process 



tr(t + 1) = T/tr(t) with rj 



1 with probability q 
-1 with probability 1 — q 
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with possible initial values tr(0) = {1,-1} both with probability ^. Notice that tr is not affected 
when (f) hits any of the two thresholds. This implies that typically new firms are created with 
positive trend. This hypothesis simplifies the analysis but does not affect the result as discussed 
in Section [51 

There are two important differences between the sign-process ^ and the trend-process ^. The 
first regards the behavior at the boundaries. Suppose both processes get to at time t — 1 
coming from a positive value at time t — 2 and remain at at time t (because, for instance, in 
the sign-process ^(i — 1) and ^(t) were negative and in the tr-process rj{t — 1) and rj{t) were 1). 
In this case, the term sign(A(/)(t)) in Eq. ([5|) is zero and therefore the sign-process will switch 
to a positive value at time t+1 with probability |. In contrast, the corresponding term tr(t) in 
Eq. ([7| can never be zero (by definition its range is {—1, +1} and the tr-process will switch to a 
positive value at time t + 1 with probability 1 — q. This means that when the noise a is small and 
therefore q is close to 1, the tr-process tends to stay longer at 0, compared to the sign-process. 
The tr-process can be easily modified to better approximate the sign-process by redefining what 
happens at zero. We discuss possibile modifications and their implications in Section [5l 

The second difference between the two processes concerns the dependencies of the draws of the 
random variables. Eq. ^ implies that sign(A(/)(t)) = sign[^(t — 1) -|- asign(A(0(t — 1)))] and 
therefore ^(t — 1) affects directly (j){t) and indirectly also 4){t + 1) through the term sign(A(?!)(t)). 
In contrast, in the tr-process the term tr(t) evolves independently of the draws of the random 
variable S, 

We now compare the tr-process with a process called persistent random walk (PRW) in the 
physics literature. PRW is a variant of the classic random walk in which the walker has a prob- 
ability q to keep the direction and 1 — q to switch direction. If we neglect the noise term a^it) 
in ^ and start with (/>(0) = 0, then (/> evolves hke a PRW on Z. The PRW obeys the telegra- 
pher's equation in the continuous limit [H [IS, [19] . An important property of the PRW is that it 
deviates, in a transient phase, from the linear scaling of the variance of the displacement with 
time, < x^ >~ t that is characteristic of the RW. Indeed, starting with all probability mass in 
zero, the variance first increases quadratically, < x^ >~ t^, due to waves that start towards — oo 
and -|-oo (ballistic scaling). After a continuous transition, the variance grows linearly as in the 
usual RW (diffusive scaling) and in the limit of large t, it evolves as j^t- Therefore, if q is close 
to 1, the variance grows still linearly for large t, although with a high diffusion coefficient j^-. 
Compared to a pure persistent random walk, our process includes, additionally, a continuous 
additive noise, a sort of reflecting lower bound at zero, an absorbing bound 6 (which leads to 
a rebirth of firms with zero fragility) , and the fact that the probability q of keeping the trend 
depends monotonously on ^. 



8M 



Jan Lorenz, Stefano Battiston: 
Systemic risk in a network fragility model analyzed with 
probability density evolution of persistent random walks 



3.2 Dynamics on the probability density function of (p 

In order to estimate the probability that the fragility (/> hits the treshold 6, we want to know how 
its pdf evolves in time, and in particular to estimate its stable pdf if this exists. 

However, it is important to notice that, at any time step t, the state of the process ^ is 
determined both by the value of 4){t) and by the value of the trend tr(t) which evolves as the 
simple two-state process ([8|). 

In order to study the evolution of the pdf of (j) one has to study the evolution of the pdf of the 
whole process (I7ll8|) 

Since the trend process takes only two values, we can divide the pdf of (/>(t) into two parts, 
corresponding to negative trend (tr(i) = —1) and positive trend (tr(t) = +1). We define the two 
functions as /"^^^ : [0,6*] -^ ]R>o and Z+^j^ : [0,6*] -^ ]R>o. 

The pdf of the whole process is determined by the pair of functions {fZf),ftu)) under the 

condition that /^ /^(t) (</>') + f^(^t)(^')d(t)' = 1. 

From this pair of functions we can derive the pdf of phi as /^(t) = fT,^^ + ftu)- In other words. 

Is' fs(t)^^')'^^' represents the probability to have fragility in [cf), (p + dcl)] and at the same time 
a downward trend, tr(t) = — 1. Analogous relation holds for the positive trend. 

It is also possible to derive the pdf of tr as ftr{t) = if fsit)^ f fld)^^ which is a pair of scalar 
values specifying the probability of having negative and positive trend and which is therefore not 
really a pdf but a probability mass function defined on {—1, +1}. 

We also define 6a to be the Dirac (5-distribution with mass shifted by a (also known as 6{- — a)), 
'*' to be the convolution operator for functions (defined for two functions /ii, /12 : M ^ M as 
{hi * h2){(p) = J hi{y)h2{^ — y)dx), f^j^ to be the pdf of the noise). 

Proposition 1. Let the pdf of 4>{t) he if^/t)^ ft/f))- If the stochastic evolution of 4> evolves as 
defined in Eq. ([7]), then the pdf of (f){t + 1) is {fl(f<i)^ fl(t+i)) with 

Qt+i) = 9rMo,e] + (b^it) + z-it))5o 

fit+i) = athofi] + (b-^it) + z+it))5o. (9) 

The functions g^ , gf are defined as 

9t = ilf^it) *5^a + {l- q)f'l(t) *^a)* f< 

gt = ((1 - Q)f^(t) * ^-» + ^/it) * ^") * f^^ (10) 

and b~{t) = f^ °° g^ , b^{t) = Jq °^ gt o-f^ the probabilities to go above for (f){t) with negative 
or, respectively, positive trend, and z^{t) = J_^g^, z~^{t) = J_^gf are the probabilities to 
go below zero for (pit) with negative or, respectively, positive trend. 
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The successive steps in the computation of g^ and gf are illustrated in Figure [2] (steps 2 to 4), 
while the computation of {fIu,i\,ft/^,-^^\) is illustrated in step 5. 

Proof. First, we look at ([9]). Notice that the convolution operation is commutative and distribu- 
tive with respect to the operation of addition. Thus, the order of the computation does not 
matter. 

Adding the noise term +cr^{t) to (/)(t) in Eq. dTj, corresponds to the convolution of the pdf of 
(p{t) with the pdf of the noise /o-g. 

The term atr(i) in the same equation, implies that the part of the pdf representing the upward 
trend is shifted upwards by a and that the part representing the downward trend is shifted 
downwards by a. This is because shifting a function along the a;-axis is represented by convolution 
with a shifted delta-function. 

If the process is on a downward trend, it will keep that trend with probability q and switch with 
probability {1 — q) . The vice-versa holds for the upward trend. Thus, a q'-fraction of fT,^. will 
remain in fZ^,i), while a (1 — g)-fraction of ftf^^ will join f7,^,-^\. The vice- versa holds for /"•". 

Finally, Eq. (fTOl) ensures that all probability mass which overlaps the interval [0, 9] is distributed 
back to [0,0]. The overlapping probability mass is determined by 6~(t), 6'*"(t), 2;~(t), 2;+(t) and 
according to the boundary conditions, it is put in a (5-peek at zero, while the trend information 
gets conserved. D 

Notice that other definitions for rebirth after failure can easily be modeled by changing Jq in Eq. 
((H) to any other pdf (for example to the pdf of the uniform distribution if firms should be reborn 
with random and equally distributed fragility). Further on, also other rules for changes of the 
trend can be modeled by replacing {b~{t) + z~{t)) and (6"'"(t) -|- z~^{t)) by other combinations. 

To better approximate the sign-process, one should replace z~{t) and z~^{t) by ^{z~{t) -|-z+(f)). 
This models the fact that a firm with fragility zero for two time steps has a zero trend, and 
switches with equal probability to the upward or downward trend, regardless of the former 
trend. 

Given an initial pdf (/7/qn, fttQ)), Proposition[T]defines a time-discrete evolution of the probability 
density function of the firm's fragility. 

In the following of this section, we will use the dynamics as defined in ^. 

Proposition 2. Consider the process defined in Eq. ^, where ^ is a normally distributed random 
variable with mean zero, variance one and pdf f^, with noise level a > 0, trend strength a > 0, 
failing threshold 9. 

If q{a,a) = Pr((T^ < a) < 1, then there exists a unique stable pdf {f^ ,f^). 
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Figure 2: Shifting, trend switching, shock addition and redistribution of overlapping mass for 
given /(^(f) = f~,^. + ftuy Parameters used: 9 = I, a = a = 0.15, the noise pdf /^^ is gaussian. 
This choice imphes q « 0.8413. In step 5 c" = (6"(t) + z-{t)) , c+ = {b+{t) + z+{t)). 

Furthermore, any initial pdf {f~,r.^, ft, r.^) converges, under the evolution defined in Proposition 
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m to (/^ ,/+) geometrically fast, with J f^ = f f^ 
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Proof. We want to apply a theorem known as Birkhoff-Jentzsch Theorem [8, Page 224, Theorem 
3]. It is an extension of the famous Perron- Probenius Theorem for nonnegative matrices to infinite- 
dimensional vector spaces. 

It is easy to see that, for any bounded pdf {fl(f.\, ft(t)) ^^^ ^^o parts of the pdf {fZfii), /J(t . lO 
are continuous on ]0,9], have a (5-peak at zero and full support [0,9]. So, after one iteration the 
dynamics ([9]) remain in the space of pairs of bounded continuous functions with a (^-peaks at 
zero. 

Let us define the operator P on the vector space of these functions such that it transforms 
{fMt)jft(t)) into {fl(^,i\,fl(t+i))- This operator fulfills the conditions of the Birkoff-Jentsch 
Theorem: it is in fact a uniformly positively bounded linear operator. 

It is bounded because, trivially, the integral of the pdf is always one. The linearity is also easily 
checked since all entities in the definition of the dynamics Eqs. (f MTU]l are linear. 
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Now we show that it is also uniformly positive (as defined in [U Page 219]). In our case an 
eigenvalue of P must be A = 1. As lower bound for (/7/^ , j^-j, /i^ , j^O we take (e,e) with e = 
(1 - g)ci(l[o,e] + (^o) with ci = f„i^{e + a). This is obviously the lowest value (/^(j+i), /^(j+i)) 
can take after one iteration because of convolution with /o-^. (Take e.g. {fl(t)ifl(t)) = (0)<^e) as 
a 'worst case'.) Further on, an upper bound exists C2(l[o,6»] + ^o) with C2 = /o-^(O). Thus, there 
exists the desired strech parameter K = C2 for the Birkoff-Jentsch Theorem. 

The Birkoff-Jentsch Theorem now states that there is a unique {f^ , f^) and that for any inital 
pdf convergence to (/~, /+) happens by iteration of the operator P geometrically fast. 

The equations j f^ = j ft = \ are obvious, because any other distribution of mass in the parts 
of the pdf would not stay constant because of the equal exchange of (1 — g) fractions in each 
step. D 

This is probably not the most general form of the theorem. Other forms of /^^ than normal (even 
with bounded support) also often lead to stabilization. But a proof is not that straight forward. 

If we exchange the terms z~{t) and z'^{t) by \{z'~ {t) + z~^{t)) to better approximate the sign- 
process, then fractions of mass in the parts of the still existing unique stable pdf will not be 
equal anymore. 

From this section we conclude that there is a unique attractive stable distribution for the prob- 
ability density of fragility in the tr process of Eq. [71 Moreover, the probability to fail at time t 

b{t) = b-it) + b+it) (11) 

converges to fixed value which we define as the limit failure probability. 

b* = lim b{t). (12) 

t— >oo 

4 Numerical results 

Unfortunately, the unique stable pdf (/+, /~) seems not to have a closed form, or at least not an 
easy one. Therefore, we compute it numerically. We set ^ = 1 (without loss of generality) and /^ 
to be Gaussian (with mean zero and variance one) and we explore the (a, o")-parameter space. 
Each pair of values {a, a) corresponds to a value of q which lies in the interval [0.5,1]. Notice 
that, assuming a different pdf for the noise would imply different values of q (cf. Section [5]). 

Figure [3] shows the first time steps of the pdf evolution for different (a, a) values. Here the initial 
value of fragility is zero and the initial value of the trend is ±1 with equal probability. Therefore, 
the initial pdf is (/ig)' fsioO ~ ^(^0,^0)- The parameter choice in the first row of plots in Figure 
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Figure 3: The first four time steps with initial pdf /^(o) = fliQ) +/(i(o) ~ f^oi-foo 
parameters, /^ is Gaussian, the g-values are computed from a and a. 



[3] corresponds to a relatively low trend strength a compared to the noise level a and thus to a 
value of q only slightly above its minimum 0.5. The random term aS, plays the major role in the 
process and in this regime the persistent random walk behaves similar to the usual random walk. 
This leads to a fast convergence of the pdf: after only four time steps (last plot in the row), the 
pdf is already close the stable pdf (cf. Figure HJ. Notice that there is a significant delta peak at 
(going beyond limit of the ordinate axis in the plot) which collects the probability to go below 
and the probability to go above 1. 
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In the second row of plots in Figure [3j the values of (a, a) correspond to values of q closer to 
one. This implies that most of the mass of the probability density function corresponding to the 
downward trend (/~) stays close to zero. On the other hand, the mass in /"•" moves with a wave 
towards the failure threshold (which is at 1, since the abscissa represents cp and 9 = 1). The 
wave smoothes out due to the repeated convolution with fa-^. Finally, in the third row of plots 
in Figure [3] g is very close to one. In this case the wave towards the failure threshold repeats 
several times until it smoothes out. Notice that in the limit a ^ 0, and thus q ^ I (not shown 
in the figure), the pdf of (j) will not converge. There will be a delta peak which moves constantly 
upwards (modulo the redistribution of its mass in zero). 

Figure H] shows instead the stable pdfs for some specific values of a and a. The pdf 's were 
computed by iteration of Eq. ^ with initial uniform distribution on [0, 9] and discretization of 
the interval [0, 1] in steps of 0.01. We proceeded until the norm of the difference in one time step 
was smaller than an accuracy level of 10~^. There were no hints that a finer discretization would 
improve the result. 

The figure shows that the stable pdf is approximatively linearly decreasing for high values of 
fragility (except for the wavy pdfs obtained with high a and low a). The slope of the linear 
decrease is non-monotonously controlled by a and a. It is easy to explain the slope in some 
cases, although this is not the case in general. When q is close to 1, it is very unlikely that a 
trajectory of the process switches direction. A trajectories with positive trend moves steadily 
along the whole range of values [0,1], repetitively hits the threshold 1 and gets reset to 0. In 
contrast a trajectory with negative trend reaches and stays there. As a result, /+ tends to a 
uniform distribution in [0, 1] and /~ tends to a delta peak in 0. On the other hand, q close to 
0.5 is implied by a much larger than a. In this regime, (j) diffuses very fast which leads again 
to a rather flat distribution for both /~ and /"*". In contrast, for intermediate values of q (for 
instance a = 0.1, a = 0.2), the profile has a pronounced negative slope for high (j). 

In the regime of high a and a close to 0, the trajectory evolves by almost discrete jumps of 
magnitude close to a. This results in a wavy stable pdf with peaks at multiples of a. But the 
wavy pattern oscillates around a line with flat slope, which is consistent with what found in the 
case of high q and sigma not too close to 0. 

We are most interested in the limit failure probability which is our proxy for the systemic risk. It 
depends on trend strength a and noise level a. So, we computed h* = 6* (a, a) for the parameter 

set a,cr g]0, 0.5]. 

Figure [6] shows that for fixed trend strength a there is an intermediate optimal a which leads 
to minimal systemic risk. In contrast, for a fixed noise level a there is no intermediate minimum 
when varying the trend strength a. Raising the trend strength always increases the systemic 
risk. The lines for high a and low a intersect. This resembles the existence of the intermediate 
optimum for fixed a. The left plot in Figure [6] shows also values of the probability to keep the 
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Figure 4: Stable pdfs for selected trend strength a and noise level a values. 

trend q at the intermediate minima of the limit failure probability with respect to o", given a fixed 
trend strength a. It turns out that the optimal noise level lies at a value of q roughly between 
0.75 and 0.9. The value of q corresponding to the local minimum decreases slowly with a. This 
is better visible in Figure [6] where we take a bird eye's view on the (a, (T)-plane, where the level 
lines of equal q appear as rays from the origin. The ordinate represents q = 0.5, the abscissa 
q=l. 

5 Robustness of results 

We checked other pdfs for the noise besides the Gaussian and in most cases we also observe 
convergence to a unique stable pdf. Notice that convergence is not assured in general by Propo- 
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a" 




Figure 5: The limit failure probability with respect to trend strength a and standard deviation 
of shocks o. (Noncomplete lines are due to extremely long convergence times.) 

sition [H We observed quantitative changes in the results but not qualitative ones in the sense 
that there always exists an optimal noise level for a fixed trend strength. 

In our model, firms fail when their fragility hits a threshold and are recreated with an initial 
value of fragility zero and an initial trend proportional to the number of failing firms with that 
trend (so mostly with upward trend) . This is a strong assumption and therefore we checked three 
other scenarios, in particular to test whether the phenomena of an intermediate optimal noise 
level is robust against these modifications. 

• If a new born firm is assigned a positive or negative trend with equal probability (instead 
of proportional to the number of failing firms of that trend) then the probability to have 
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b (a,a) 




Figure 6: Contour plot for the limit failure probability s{a,a) and contour lines for q{a,a) = 
0.8, 0.9. The solid black line denotes the optimal value of a regarding a fixed a. 

a positive trend / ff,^<. converges to a fixed number below | which depends on q. In the 
extreme case, q = 1 it goes to zero. That would implies that the probability to fail will also 
go to zero in the limit. We saw that for a fixed trend strength there is a critical noise level 
that implies such high q that the systemic risk drops to zero when the noise level gets below. 
Nevertheless, for low trend strength values (a below about 0.12) an intermediate optimal 
noise level still exists until further decreasing the noise level causes the sudden drop due to 
the extinction of the upward trend. One may criticize this variation of the model because 
it does not converge to equal proportions of positive and negative trend. But stable equal 
probabilities for upward and downward trend seems quite reasonable because judgement of 
fitness is always done comparatively in an economy. If economy divides firms in good and 
bad ones this should not lead to a possible die out of one class. 

• Another suggestion against our original model could be that firms are not born with zero 



IT/EO] 



Jan Lorenz, Stefano Battiston: 
Systemic risk in a network fragility model analyzed with 
probability density evolution of persistent random walks 



fragility but i.e. random and uniformly distributed in the fragility interval. This obviously 
changes the limit pdf, but at least in this example the qualitative behavior with the exis- 
tence of an optimal number of hedging partners for given rend strength remains the same. 

• Another idea is to renormalize the probability mass after a failure. We do this as follows: 
we do not redistribute the probability mass after a failure to zero but just rescale /"*" 
proportional to its actual shape such that it has the same total amount as before. The 
same with /~. On the level of individual firms this means that new firms are born with 
fragilities drawn randomly from the actual distribution of fragilities with that trend. That 
means if the distribution of fragilities is double peaked new firms are most likely to appear 
with fragilities around that two peaks. This dynamics imply that a given peak structure 
gets amplified by the evolution of new firms. In fact this dynamic fragility distribution for 
new firms leads to an amplification of mass in high fragility intervals. That means that 
with high probability new firms are born with high fragility (which seems reasoable). In 
the limit these regimes are characterized by virtually all firms with positve trend failing 
each year. That means that decreasing the noise level (which increases q) is even more 
dangerous. Nevertheless, there still exists an intermediate optimal noise level for a given 
trend strength to minimize the systemic risk. 

6 Conclusions 

We have presented a simple model for the stochastic evolution of the fragility of units in a network. 
The model applies in particular to networks of firms connected via financial relationships. The 
basic ingredients of the model consist in a mechanism of risk sharing that leads to decrease 
the fluctuation of the fragility and in a mechanism of reinforcing feedback on the fragility from 
the trend in the immediate past of the fragility of the firm itself and its neighbors. Under this 
assumptions, the number of bankruptcies in the system is minimized for an intermediate density 
of links in the network e.g. for an intermediate number of hedging partners. The result is of 
interest from the point of view of policy design for the control of systemic risk. 

The effect depends strongly on a dynamics divisions of firms into two classes: the good evolving 
(with decreasing fragility) and the bad evolving firms (with increasing fragility). One might 
question that this hard cut between the two classes exists. But we argue that actually, slight 
differences in performance are exacty what investors like hedge funds search for when they try 
to profit from investments indepently of the economic trend. So, even very slight differences may 
matter a lot for reinforcing trends. Further on, these kind of investment strategies have become 
more popular. 
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With respect to the original model, the analysis presented here neglects the process of cascades of 
failures and therefore underestimates the number of joint failures. However, its advantage is that 
the evolution of the probability distribution of failures can be expressed analitically and that the 
stable distribution (which we prove to exist and be unique) can be computed numerically. 

The impact of heterogeneity in the topology of the network is not studied at this stage. Further- 
more, the hedging network is not dynamic. This implies for instance that firms do not have the 
possibility to interrupt hedging relations with partner who do not perform well. This assumption 
is certainly not very realistic on a time scale of years. However, it is also true that many partner- 
ship or insurance contracts cannot be modified in a very short time. Furthermore, in future work 
the impact of heterogeneous trend strength, noise level and failing threshold should be studied. 
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